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Abstract 

We present numerical results for the statistics of z's (a's are defined as 
logarithm of eigenvalues of the transfermatrix T^T) at the critical points 
of Anderson transition in 3D and 4D. The change of the density of z due to 
the crossover from the metallic to the localized regime is described. Linear 
behavior p{^z) — z sX the critical point in 3D is proven and discussed. In 
the insulating regime, the universal form of p has been found. 

While the finite-size scaling analysis of the disorder-induced metal-insulator 
transition enables us to find critical disorder Wc and critical exponent v from 
the knowledge of the smallest positive Lyapunov exponent, z\ in the quasi-one 
dimensional limit |l| , the knowledge of the higher z's are necessary to understand 
the statistical properties of transport in cubic samples. Here, Zi is defined as 
logarithm of the i-th eigenvalue of the matrix T^T (T is the transfer matrix). 

Statistics of z's is well known understood in the limit of the small disorder. 
Its properties, derived in the analogy with the random matrix theory, enables 
us to explain universal features of the transport in weakly disordered mesoscopic 
systems The key role in this explanation plays the density defined as 

p(z) = (^<5(z-zO), (1) 

where the summation covers all channels and brackets mean averaging over 
statistical ensemble. In the weak disorder limit, p(z) = const. Another typical 
characteristics of the distribution [5) of (normalized) differences between z's. In 
the weak disorder fimit, P{S) equals to the Wigner surmise P(s) = exp — -Is^ 

i- 

Our aim in this paper is to shown how the statistics of z's changes when 
system undergoes the metal-insulator transition. We believe that the under- 
standing of the statistical properties of z's at the critical region would provide 
us with serious basis for the more general understanding of the Anderson tran- 
sition, including the description of the system size and disorder dependence of 
conductance and its statistical moments in the critical region. 

We consider Anderson model: 

7^ = iy^e„|n)H+ ^ |n)(n'|. (2) 

n [n.n.] 

In (|^), 71 numerates sites on the c?-dimensional cubic lattice L'', and W measures 
the disorder. For the box distributed random energies e, |e| < 1/2, model 
exhibits the metal-insulator transition at the critical point Wc ~ 16.5 (34.5) for 
d = 3 (4) , respectively. 
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Figure 1: P{5) for 3D (full symbols) and 4D (open symbols) at the critical point. 
Solid lines are Wigner surmise and Semipoisson distribution. 

Recently, [Q, |^ some speculative models has been proposed, in which weak 
disorder statistics of z's has been generalized to models which reflect some spe- 
cial features of the metal-insulator transition. Unfortunately, no one of them 
succeeded to describe the metal-insulator transition completely. 

The analysis of the statistics of z's has its counterpart in the level statistics 
. Here, the most important quantity is the distribution P{s) of the (normal- 
ized) differences between energy levels. It has been shown, that there are three 
typical form of P{s): Wigner surmises (WS) for the metallic. Poison for the 
insulating and the third, universal distribution at the critical point. The same 
scenario has been found for the statistics P{S) of the (normalized) differences 
between z's 1^. 

We present in Fig. 1 P{5) in the critical point. Data show that (i) P{6) is sys- 
tem size invariant, (ii) it depends on the dimension, and that (iii) it follows nei- 
ther the Wigner surmises nor the Semipoison distribution P{s) — 4sexp(— 2s). 

The first too results correspond to that of the level statistics; the second one 
confirms that the statistics of z's is more similar to the localized one in higher 
dimension However, the third one is in the contradiction with numerical 
observation reported in Ref. |^ , where the Semipoisson distribution of the level 
statistics has been found as the result of the sensitivity of the level statistics 
to the boundary conditions at the critical point. Although we consider only 
periodic boundary condition in this work, we do not believe that the use of dif- 
ferent boundaries will remove this disagreement. The reason is that the typical 
values of z's are rather large at the critical point: The distribution of higher z's 
is Gaussian with mean value > 5 and variance ~ 1 It is therefore highly 
unprobale that the change of boundaries could infiuence the statistics of higher 



3 



8 



3D Anderson model 
W=6,1 0,1 6.5,25,32 




T 



6 



T 



P(z) 




60 



Figure 2: The density p{z) for cubic system with L = 12. 



z. In other words, higher channels are strongly localized (although mutually 
correlated) and so they do not feel the form of the boundaries. 

Fig. 2. shows how the density p{z) changes due to the growth of disorder. 
The typical semicircle form of the p{z) could be found only for extremely small 
disorder {W = 6) (we remove the contribution of closed channels in this case). 
The growth of the disorder deformates p{z): the last decreases for small z in 
favor of the maximum, which moves towards the higher values of z. At the 
critical point, p{z) becomes linear (quadratic) in z in 3D (4D), respectively. 
This agrees with our previous result, obtained for the quasi-one dimensional 
samples 

The critical density deserves the special Figure. In Fig. 3a we present p(z) 
for different size of the sample. It is evident, that only small part of p(z) is 
system size invariant. The interval of the size invariance equals approximately 
to that at which p{z) remains linear. When comparing p{z) for different system 
sizes, > iminj then the interval, at which all p's coincides, is z < Lniin/2. This 
restriction must be taken into account in analysis of the scaling properties of 
higher z's 0. 

In Fig. 3b, c we present the normalized critical density for 3D and 4D An- 
derson models. For 4D, p{z) ~ in the limit of small z. 

The form of p{z) remain approximately the same also when W exceeds crit- 
ical point. The whole distribution only shifts towards the higher values of z. In 
the limit of large disorder we find that 



(Fig. 4). The most remarkable differences could be found only in the tail of the 
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Figure 3: (a) Critical p{z) for different system size. Only small (linear) part of 
p{z) is system-size invariant, (b) normalized p{z) for 3D and (c) 4D system. 

density for small values of z. 

In conclusion, we present numerical data for Lyapunov exponents for 3D and 
4D Anderson models in the neighbor of the critical point. We show how two the 
most famous characteristics of the statistics of zs' change when the strength of 
the disorder exceeds the weak disorder limit. In difference to the level statistics, 
we found no Scmipoisson distribution of the differences of z's. Our data support 
belief that the statistics of z's and, consequently, of the conductance has a simple 
form similar that developed for the metallic regime. 
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Figure 4: Density p{z) for the 3D system with W = 32 (o), 45 {<}) and 60 
(fuU symbols) for two system size {L = 8 and 12). We shifted the densities for 
W ^32 and 45 in {zi{W = 60)) - {zi{W)) to show their the universal form (|). 
(zi) for = 60 is 20 (3) for L = 8 (12), respectively. 
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